PROCEEDINGS 

OF SCIENCE 



in 
o 
o 

(N 

> 
O 

o 



> 
00 
OS 

o 
in 
o 



X 



Parameter dependence of the topology change and 
the scaling properties of the topology conserving 
gauge action 

KEK-CP-169 

YITP-05-54 

KUNS-1987 

Hidenori Fukaya, Tetsuya Onogi* 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan. 

Shoji Hashimoto 

High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801, Japan. 
School of High Energy Accelerator Science, The Graduate University for Advanced Studies 
(Sokendai), Tsukuba 305-0801, Japan. 

Takuya Hirohashi 

Department of Physics, Kyoto University, Kyoto 606-8502, Japan. 

Kenji Ogawa 

School of High Energy Accelerator Science, The Graduate University for Advanced Studies 
(Sokendai), Tsukuba 305-0801, Japan. 

The topology conserving gauge action proposed by Liischer is expected to reduce the number of 
non-smooth gauge configurations as well as the topology change compared to the conventional 
actions for the same lattice spacings. We report our quenched QCD study of the topological 
stability and the scaling violation of the static quark potential. We find that the the topology 
change is indeed suppressed when the parameter e is of order one. We also find that the scaling 
violation in the static quark potential remain reasonably small in the parameter range of our study. 
Our study is done at the inverse lattice spacing a = 1.4-2.5 GeV with the lattice size L = 1.0- 
1.6 fm. 
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1. Introduction 

The admissibility condition of the gauge fields on the lattice ||] 

||1 — //i V (jc)|| < £ for all x, jj,, v, (1.1) 
guarantees the locality and the uniqueness of the index of the overlap-Dirac operator. [0, f|]. As an 



example, the gauge action which restricts the gauge field to satisfy the bound (1.1) 



_ l-ReTr/y(x)/3 
P ^l-(l-ReTrP MV (x)/3)/e' U ^ 

was proposed by Luescher [JTJj . This action may be useful for QCD simulations in two reasons: (1) 
it may serve to efficiently collect gauge configurations with fixed topology Q ^ in the e-regime, 
(2) it may give a possibility to reduce the numerical cost of dynamical overlap fermions with fewer 
low-lying modes of Hw and less frequent topology change. In fact, this action has been proven to 
be useful for the massive Schwinger model for stabilizing the topological charge and for improving 
the chiral behavior of the domain wall fermions [g, gj. Also in the four-dimensional quenched 
QCD, the good stability of the topological charge has been observed with some reasonable choices 



of parameters |0 g|. See also [|10|]. 

In this report we present our quenched study of the topology conserving gauge action ( |1.2[ ) on 
the stability of the topology and the scaling. Studies on the low-mode distribution in quenched and 



dynamical QCD are reported by Matsufuru [11]. 



2. Lattice setup 

In our study, we take three values of l/e (= 0, 2/3, and 1) with three lattice sizes 12 4 , 16 4 , and 
20 4 and four lattice spacings in the range aT x = 1.4-2.5 GeV. Here, l/e = corresponds to the 
conventional plaquette action and l/e = 2/3 is the boundary, below which, all gauge configurations 
are allowed when the gauge group is SU (3). Link variables are generated by the standard hybrid 
Monte Carlo algorithm with At = 0.01-0.02 and N m d s = 20^-0, where N mc i s is the number of 
molecular dynamics steps and At denotes its step-size. We accumulated at least 2,000 trajectories 
for thermalization before measuring observables. We monitored the plaquette values of the gauge 
fields, but we did observe any case where the admissibility condition is violated through the hybrid 
Monte Carlo updates. 



3. Stability of the topological charge 

In order to measure the topological charge, we developed a new cooling method, in which we 
carry out the hybrid Monte Carlo simulation using the topology conserving action with an expo- 
nentially increasing coupling j8 = j3 coo i an d an exponentially decreasing step size At as functions of 
trajectory n t . After 50-200 steps, the gauge fields are expected to be cooled down close to the clas- 
sical background in each topological sector. In fact, the geometrical definition of the topological 
charge [ |I2| ] 

Qseo = ^"-2 £^ vpff ReTr(P MV (x)P pff (x)) (3.1) 
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Figure 1: Stability of the topological charge (3.1) for each parameter set. Results are plotted as a function 
of lattice spacing squared. Three values of 1 /e are distinguished by lines: (dotted), 2/3 (dashed), 1 (solid). 
Circles, squares, and up-triangles correspond to lattice sizes 12, 16 and 20, respectively. The down-triangles 
are from M\ . 



of these "cooled" configurations gives numbers close to an integer. 
Figure [l] shows the stability of the topological charge defined as 



Stab G = 



trj 



Tplaq x $Q 



(3-2) 



where T p i aq is the autocorrelation time of the plaquette, measured by the method in Appendix E 
of JT3|], Af tr j denotes the number of trajectories and #Q is the number of topology changes. Since 
the topological charge is measured every 10-20 trajectories, Stabg gives only an upper limit. Our 
results show that the topological charge becomes more stable for higher l/e, smaller L and smaller 



4. Static quark potential 

We calculate the static quark potential V (r) from the Wilson loops W(r,t) for every 20 trajec- 
tories. The spatial side r is taken to be in 6 different directions of the 3-dimensional unit vectors: u 
= (1,0,0), (1,1,0), (2,1,0), (1,1,1), (2,1,1), (2,2,1). We calculate the Sommer scales r and r c defined 
as r^F{ro) = 1.65 and r^F(r c ) = 0.65, respectively [[14], 15]. F is the force obtained from the static 
quark potential. We choose a Q = configuration for the initial condition of HMC steps and do not 
take care of the topological charge hysteresis assuming that the topology of the gauge fields would 
not affect the Wilson loops if they are small enough. 

Figure § shows the scaling of r c /ro. One can see that they agree well with the results with the 
plaquette action and its continuum limit JT5| ] in the region (a/ro) 2 < 0.08. Also, as seen in Figure || 
for long distances, the quark potential itself V(r) = ro(V(r) — V(r c )) does show a good agreement 



with that in the continuum limit obtained with the plaquette action [15]. Here, V{r c ) is measured 
using an interpolation polynomial of order 5. For short distances, they show 10-20% deviations 
due to the violation of the rotational symmetry but the discrepancies are comparable to those with 
the plaquette action. 
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Figure 2: Scaling of r c /ro. Squares and triangles are data for 1/e = 2/3 and for 1, respectively. Open 
circles are the result with the plaquette action in jl5[], while the filled circle denotes their continuum limit. 
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Figure 3: Left: heavy quark potential at j3 = 1.3, 1/e = 1.0 on a 12 4 lattice. Dashed line shows the 
continuum limit obtained by an interpolation from the result of jl5|]. Different symbols show the V(r)'s 
with different orientations parallel to w's. Right: (V(r) -~Vcont{r))/V C ont{r), where V C ont(r) represents the 
continuum limit. The error of V cont (r) is ignored (< 1%). 



We conclude that the heavy quark potential for the admissible gauge fields is quite similar with 
that of the plaquette action and no serious inconsistency can be found. 



5. Perturbative renormalization of the coupling 

Ellis and Martinelli computed the two-loop corrections to the gauge coupling for general gauge 
actions which can be written by the plaquette Jl6|]. Using their formula, the renormalized gauge 
couping gM in the so-called Manton scheme is expressed by the bare coupling go as 



1 



«Af(l/«) #0 



+M+M80, 



(5.1) 



where (A 1} A 2 ) = (-0.2083,-0.03056), (0.34722,-0.04783), and (0.625,-0.10276) for 1/e = 0, 
2/3, and 1 , respectively. In Figure ^ we plot the inverse squared coupling in the Manton scheme at 
a reference scale jj. = 5/ro as a function of the lattice spacing. Here, we use the one-loop formula 
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Figure 4: 1 /glf(5/ro) in the Manton scheme as a function of lattice spacing squared. Results for various 



1/e and r$ are plotted. The bare coupling 1 /g^ is also shown by open symbols for comparison. 



in ( p.rj ) for the change of scheme and two-loop renormalization group equation for the evolution 
to the reference scale. We find that the renormalized coupling for various values of 1 /e shows a 
good agreement, which means that the change of 1 /e is well described by the perturbation theory. 
One can also see that the scaling violation is small for the renormalized coupling. However, we 
note that the two-loop result gives larger corrections for larger 1 /e so that the convergence of the 
perturbative series are poor for 1/e = 1 and marginal for 1/e = 2/3. 

6. Summary 

We studied the topology conserving gauge action in the quenched approximation for various 
values of j8, 1/e and the lattice size. By measuring the topological charge with a new cooling 
method, we find that the stability of the topological charge is improved for larger 1/e when com- 
pared at a same lattice spacing and lattice size. Measuring also the static quark potential and the 
Sommer scales r^, r c , the scaling violation are found to be small. Our study shows that the topology 
conserving gauge action is feasible for QCD simulations. Practical applications of this action such 
as QCD in e-regime and simulations with dynamical overlap fermions are underway. 
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